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23 ■ ABSTRACT 



The validity of a photon shower generator QEDPS has been examined in detail. 
This is formulated based on the leading-logarithmic renormalization equation for 
the electron structure function and it provides a photon shower along the initial 
e^. The main interest in the present work is to test the reliability of the generator 
to describe a process accompanying hard photons which are detected. For this 
purpose, by taking the HZ production as the basic reaction, the total cross section 
and some distributions of the hard photons are compared between two cases that 
these photons come from either those generated by QEDPS or the hard process 
e + e~ — > HZ 0/ y ■ ■ • 7. The comparison performed for the single and the double hard 
photon has shown a satisfactory agreement which demonstrated that the model is 
self-consistent. 
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1 Introduction 



In a series of works |l|]-|§ we have proposed a generator QEDPS which develops 
a photon shower along the initial e of any hard annihilation process. This gives 
the correction by the initial state radiation(ISR). Usually ISR is implemented by 
the structure function of the initial e . By using the generator, however, it is 
possible not only to give the corrected total cross section with arbitrary cuts but 
also a distribution over any kinematical variable including transverse momentum of 
photons. It should be also noticed that the annihilating e + and e~ no more make 
the head-on collision, as they deviate from the beam axis by the radiation. 

In the following sections we examine in detail the validity of QEDPS. We take 
e + e~ — > HZ together with possible hard photons as the hard process. The next 
section will make a brief review of QEDPS. We also summarize the analytic form of 
the structure function up to 0(a 2 ), which should be compared against the photon 
shower in section 3. Section 4 is devoted to the numerical test of QEDPS. The 
correction to the total cross section is tested in the sub section 4.1 between the 
simulation and the analytic structure function. The subsequent two sub sections, 
4.2 and 4.3, are devoted to a single and double hard photon emission, respectively. 
The last section gives the summary and discussions. 



We review briefly the formulation of QEDPS in this section. The basic assumption is 
that the structure function of an electron, with the virtuality Q 2 and the momentum 
fraction x, obeys the Altarelli-Parisi equation 



2 QEDPS 
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Here e is a small quantity specified later. Then the original equation can be converted 
to the integral equation. 

a rQ 2 dK 2 r 1 -* dn 

D(x, Q 2 ) = H(Q 2 , Ql)D(x, + — — -H(Q 2 , K 2 ) / ^-P(y)D(x/y, K 2 ). 

iix jQi K z Jx y 

(3) 

Here, rigorously speaking, Q 2 should be m 2 e as it gives the initial condition. For sim- 
plicity the fine structure constant a is assumed not running with Q 2 . The function 
n, which is nothing but the Sudakov factor, is given by 

nw^Hex P (-^/;:^/;-w w ), ( 4) 

and implies the probability that the electron evolves from Q' 2 to Q 2 without emitting 
hard photon. In other words n already contains the contributions from the soft 
photon emission, which causes the change of the electron virtuality, and from the 
loop corrections in all orders of perturbation. 

The integral form Eq.(|3|) can be solved by iterating the right-hand side in a 
successive way. Then it is apparent that the emission of n photons corresponds to 
the n-th iteration. Hence it is possible to regard the process as a stochastic one that 
suggests the following algorithm of the photon shower]]]]. 

(a) Set Xb — 1. The variable Xb is the fraction of the light-cone momentum of the 
virtual electron that annihilates. 

(b) Choose a random number t]. If it is smaller than H(Q 2 , Q 2 ), then the evolution 
stops. If not, one finds the virtuality K 2 that satisfies rj = Ii{K 2 , Q 2 ) with which a 
branching takes place. 

(c) Fix x according to the probability P(x) between and 1 — e. Then Xb is 
replaced by XbX. One should go to (b) by substituting K 2 into Q 2 and repeat until 
it stops. 

Once an exclusive process is fixed by this algorithm, each branching of a photon 
in the process is dealt with as a true process, that is, an electron with x, K 2 decays 
as 

e-(x, -K 2 ) - e-(xy, -K' 2 ) + 7 (x(l - y), Q 2 ). (5) 
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In the infinite momentum frame where the parent has the momentum p = (E, Ot, p z ), 
E = yjp 2 z — K 2 , p z = xp*, the daughter electron and the photon have p' = (E 1 , k T , yp z ) 
and k = (k , — k r , (1 — y)p z ), respectively, with E' = \J y 2 p 2 z + — K' 2 and k = 
\f (I — y) 2 P 2 z + + Qo- Here we have introduced a cutoff Ql to avoid the infrared 
divergence. The momentum conservation at the branching gives 

- K 2 = -K' 2 /y + Ql/{1 -y) + 4/(y(l - y)), (6) 

which in turn determines k^ from y, K 2 , K' 2 . Hence one can get the information 
on the transverse momentum. This enables us to give the k 2 -. distribution by the 
simulation as well as the shape of x. 

Further the kinematical boundary y{K 2 +Q"^/ [l — y)) < K' 2 , equivalent to ky > 0, 
fixes e as 

e = Ql/K' 2 , (7) 

since K 2 K' 2 is expected. In ref. it has been demonstrated that the existence of 
e which depends on the virtuality gives rise to a small but non-negligible contribution 
to the structure function, though it is a constant. This is because when the virtuality 
becomes to the order of Qo, e is not small, while in Eq.@ it is assumed to be small 
enough. 

The above description of the algorithm has shown the single cascade scheme. This 
implies that either of e~ or e + is able to radiate photons when the axial gauge vector 
is chosen along the momentum of the other electron, namely e + or e~. In writing a 
computer code for the shower, however, it is convenient to use the double cascade 
scheme to ensure the symmetry of the radiation between e + and e~[Q. It can be 
shown that these two are mathematically equivalent in the LL approximation. The 
extra finite contribution due to e is, however, different between single and double 
cascade scheme. The details will be found in ref. 0. 

In the formulation there are two parameters Q 2 and Qq. In the program the 
following values are chosen: 

Q 2 = m 2 e e = m 2 e x 2.71828 Q 2 = 1(T 12 GeV 2 . (8) 
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The former value was settled to take into account effectively the constant term —1 of 
(5 in such a way j3 = (2a/7r)(hi(s/m 2 ,) — 1) = (2a /tt) \yi(s / (m 2 e e)) . Since the second 
parameter is unphysical, any physical observable should not depend on it. It has 
been checked that increasing Qq up to O(ml/10) leaves the result unchanged in the 
statistical error of the event generation. 



3 Structure function 



The differential equation Eq.(TI|) can be easily solved if one introduces the moment 
of the structure function byH 



D(n,s) = f 1 dxx n ~ 1 D(x,s). 
Jo 

The solution is then 

D(n, s) = exp 
with the initial condition D(x, Q 2 S ) = 5(1 — x). Here 



2 V2 u J n n + 1 



Si(n) = £ 



(9) 



(10) 



is an analytic continuation of the finite sum Dj=i 1/0 into the complex n-plane. We 
expand Eq.([U]) with respect to (3 in the following way 



D(n,s) 



1 + ^ + ^ 2 ) exp (-f s ' ( - i; 



4a + 



4 / 2 \n n + 1 



4 J 2 \n n + 1 
and further for n — > oo 



(12) 



exp f-^Si(n - r 



exp (—(ie + Inn) j - ^[Si(n - 1) - 7s - Inn] 



■5 



+ 



Iff 

2! I 2 



{[Sl(n _l)]2_ (7s + lnn) 2 }; 



(13) 



where 7^ = 0.57721 • • • is the Euler's constant. The inverse transformation defined 
by 



D(x,s) 



c + io ° dn 



x- n D(n,s), 



1 c—100 2,7li 

with c > being a real number to fix the integration path, gives 



(14) 



D(x,s) 



9 

128 

2 

'-(1 + x) - — 
A y ! 32 



1 + / + 



C(2)' 



/? 2 



2 V 



,,8/2-1 



4(1 + z) ln(l-x) + 



1 + 3x 2 
1 — x 



lnx + (5 + x) 



(15) 



Here we have made a resummation to get the first term (1 — x)' 3 / 2-1 . 

Eq.flTo]) is used throughout for the test of QEDPS in the following sections. One 
should keep in mind that these results are obtained in the leading-log(LL) approxi- 
mation and non-leading terms are not fully included corresponding to QEDPS. 



4 Numerical test of QEDPS 
4.1 Corrected total cross section 

In Ref. HI we have already compared the model with the analytic structure function 
and also with the perturbative calculations of the initial state radiation(ISR) avail- 
able up to O(a 2 )0. In these comparisons we had obtained a satisfactory agreement, 
but as explained in section 1 the recovery of the missing finite term should give bet- 
ter agreement. This is important to be confirmed since the experimental accuracy 
is becoming less than 1%. 
The basic hard process is 

e + e' -> HZ . (16) 

We take this reaction as our example since it is not only the most suitable one to test 
the ISR corrections as there is no final state radiation but also very simple process 
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with a few Feynman diagrams even when some extra photons are attached. Its 
cross section a"o(s) is combined with QEDPS to generate the events with indefinite 
number of photons. None of them is assumed to be observed, that is, no cut is 
imposed on the photons. On the other hand the corrected total cross section can 
be easily estimated by convoluting <7o(s) with the structure function. In this work 
all the hard process is calculated by Monte Carlo integration 0. In other words we 
also have a generator for the hard process. 

In Table 1 we show the results evaluated by these two different methods. The 
used parameters are given by 



a = 1/137.036, m e = 0.511 x 10~ 3 GeV, 
M w = 80.230 GeV, M z = 91.188 GeV. 



(17) 



The mass of Higgs particle is taken to be the same as M^y. This has no realistic 
meaning. The total energy is chosen in the range \/s = 200 — 1000 GeV, roughly 
from LEP200 to linear colliders. From the table one finds that the simulation agrees 
well with the analytic formula within the accuracy less than 0.4%. This indicates 
that QEDPS effectively reproduces the structure function which includes the soft 
photon resummation and also the expansion of terms up to 0{f3 2 ). 



v^(GeV) 


HZ/ps 


HZ /si 


(#Z/ps)/(#Z/sf)-l 


200 


(5.727±0.004)xl0- 1 


(5.727±0.002)xl0- 1 


0.0% 


300 


(2.411±0.002)xl0- 1 


(2.419±0.001)xl0- 1 


-0.33% 


400 


(1.193±0.002)xl0- 1 


(1.196±0.001)xl0- 1 


-0.25% 


500 


(7.125±0.007)xl0" 2 


(7.137±0.003)xl0- 2 


-0.17% 



Table 1 Total cross sections(pb) without cut. 

Finally we compare QEDPS with the fixed order perturbative calculations for the 
higher order corrections. At ^/s = 500 GeV, for example, the total cross sections are 
a = 5.766 x 10~ 2 pb in the tree level (Eq. fllfiD), U\ = 7.181 x 10" 2 pb for the corrected 
one up to 0(a) and a 2 = 7.208 x 10~ 2 pb up to 0(a 2 )|| |10[, respectively. The 



numbers in Table 1 do not include the so called .fT-factor (K = l+a/7r(7r 2 /3 — 1/2) = 
1.00648...) which recovers almost the constant term in 0(a). Multiplying this value, 
we get 7.171 x 10~ 2 pb of HZ/ps or 7.183 x 10~ 2 pb of HZ /sf, thus the convergence 
of the perturbation series shows a good behavior. 



4.2 Single hard photon test 

Let us consider the following experimental situation. Together with H and Z° we 
shall observe a hard photon with E 7 > 1 GeV outside of the cone of 5° from the 
beam axis. Inside the cone any number of photons can be radiated. The photons 
with E~ < 1 GeV are regarded as soft and get no angular limitation. Hence one 
detects HZ ^ as the final state. 

To estimate the total cross section and distributions with respect to some kine- 
matical variables, which we specify later, we consider four different ways to estimate 
this cross section. The first one is simply to calculate it in the tree level 

e + e- -> HZ ^. (18) 

Obviously this gives the lowest order result without any radiative correction. It is 
not expected to be precise enough then, but we include this result in our comparison. 

The second one is the same process but corrected up to 0(a) in perturbation, 
which should give more accurate estimation. In this case one has to evaluate the 
cross section for the process 

e+e" -> #Z° 77 , (19) 

in the tree level under the experimental condition on the photons as same as before. 
Also one needs the one-loop correction to the process Eq . (|i~8l) . All the necessary 
calculations are found in Ref. |§. 

The third way is to impose the condition on the photons generated by QEDPS 
being combined with dao(s). This is the same thing as given in the previous section 
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but with the cut. We denote this estimation as HZ/^ps to imply that QEDPS is 
applied to the bare process Eq.flTBD but one of the generated photons is observed. 

The last one is to dress the photon shower to the radiative process Eq.(|T8"D in 
the tree level. In this case the photon associated with HZ production is regarded 
as the observed one. Other photons supplied by the simulation are supposed to be 
invisible so that they escape inside of the cone or are soft. We denote this process 
as HZ'y/ps. Here it is important not to make double counting of the hard photon. 
This can be achieved by the following manner. In Eq.fl^) the virtuality of the initial 
e ± is required to satisfy K 2 < K' 2 at each branching in the cascade. This order 
should be also maintained for the electron in the process Eq. (|T8|) . In other words 
the virtuality of in the cascade is restricted to be smaller than that appears in 
the hard process. 

One should notice that the comparison of the third and the fourth cases provides 
a self-consistency check of the model while the first two allows the comparison with 
the perturbative calculations. 

The results of the total cross sections are summarized in Table 2 which are ob- 
tained by these four different methods at various energies with the cut. 



v^GeV) 


HZ'y (tree) 


HZ j (up to 0(a)) 


HZ/^ps 


HZj/ps 


200 


5.602xl0~ 2 


4.418X10" 2 


(4.534±0.008)xl0-2 


(4.424±0.002)xl0-2 


300 


3.284xl0~ 2 


3.018xl0~ 2 


(3.086±0.016)xl0" 2 


(2.945±0.002)xl0~ 2 


400 


1.765xl0- 2 


1. 701xl0- 2 


(1.754±0.010)xl0-2 


(1.648±0.001)xl0-2 


500 


1.101xl0~ 2 


1.088xl0- 2 


(1.143±0.007)xl0-2 


(1.048±0.001)xlQ-2 



Table 2 Total cross sections(pb) : single-hard photon required. 



At 500 GeV, the effect of 0(a) corrections is small, less than 1%. It does not 
mean, however, that Born approximation is enough because the various distributions 
are different between them as shown in Fig.l. The plots show the first case and 
the histograms does the second one. Q is the virtuality of the s-channel virtual 



9 



boson(Z°) or the invariant mass of the HZ system. Q z is the longitudinal momenta 
of this system. E 1 and E] are the energy and the transverse momentum of the 
required hard photon with respect to the beam axis, respectively. There is seen a 
large discrepancy in the Q and Q z distributions. This is not surprising because there 
is no correction from the soft photon emission in the tree level case. On the other 
hand, E 1 and E] distributions agrees rather well. 

Next we show in Fig. 2 the comparison of the same distributions between the third 
case, HZ/^fps and the second case, HZ'-f with the 0(a) corrections. Here, again, 
the histograms show the second case and the plots the third one. This time Q and 
Q z distribution are quite consistent because the both cases include the soft photon 
emission. In E 1 and Ej distributions the third case shows an enhancement against 
to the second, but still consistent. One reason of this enhancement is presumably 
attributed to the fact that the constant terms of 0(ajn) are missing in the LL 
approximation. Because of this enhancement, the total cross section of HZ/'jps are 
always larger than that for HZj with the 0(a) corrections. QEDPS model, however, 
gives the events with the acurracy in the level of a few percent even though a hard 
photon is required. 

Fig. 3 shows the same comparison as in Figs.l and 2 but between the fourth case, 
HZ'-f /ps, given by the plots and the second, HZ'-f with 0(a) corrections, by the 
histograms. The agreement of all the distributions is very well. But the total cross 
sections in Table 2 of the fourth case are always smaller than the second one. This 
corresponds to the fact that in the lowest edge of the E} , the differential cross 
section of the fourth case is smaller than the second one. As mentioned above, we 
require the order of the virtuality to combine the one photon emission process with 
the QEDPS model. This recipe may be still naive and has some ambiguity in the 
very small El region. On the other hand, in the large El region, it works well. As 
long as one intends to study the events with a hard photon, the fourth is the best 
one. 

From these results one can conclude that QEDPS is consistent not only with the 
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result of perturbation calculation up to the 0(a) corrections but also within itself. 
This is a strong support for the reliability of QEDPS. 

What happens if one convolutes the structure function with the radiative process 
Eq.fll8l), that is, HZj/sfl It is obvious that one has the double-counting because the 
hard photon is also contained in the structure function. Hence a naive convolution 
leads to an overestimation and one should introduce some special prescription to 



get a correct result [11 



It would be, however, interesting to see if QEDPS could reproduce the naive 
HZj/sf. In the former, as explained in the previous section, the virtuality cut 
prevents the double-counting. Then one may expect that to remove this cut would 
give the same result as HZ~//sf. One kinematical difference should be noticed. In 
the case of the structure function the head-on collision of e takes place. On the 
other hand, with QEDPS the annihilation of electrons in the hard process cannot 
make the head-on collision because an emitted photon has a finite pj. To compare 
these two on the same ground then, we require the hard photon cut in the CM 
system of the annihilating e . The results summarized in Table 3 show that the 
above consideration is legitimate. The total cross section in the first row calculated 
by the structure function does not show any significant difference between Lab. and 
CM system while that by QEDPS in the second row considerably changes, resulting 
the consistent number with HZ'y/sf only where the cut is imposed in the CM system. 
These results are unphysical after all, as the double-counting is allowed and the cut 
in the CM system is artificial. The comparison, however, provides another evidence 
that QEDPS is consistent with the structure function. 



hard photon cut 


in Lab. sys. 


in CM sys. 


a(HZj/sf) 


(1.256±0.002)xl(r 2 


(1.266±0.001)xl(T 2 


a(HZ^ /ps) without 
the virtuality cut 


(1.415±0.005)xl(T 2 


(1.260±0.001)xl(r 2 



Table 3 Total cross sections(pb) at y / i=500GeV : single-hard photon required. 



11 



4.3 Double hard photon test 

Further check of the model is possible by requiring the double hard photon emis- 
sion. In this case, however, the perturbative calculation is restricted to only the 
tree process, Eq.fll^), as we have no loop-correction of 0(a 3 ) yet. Hence the first 
estimation comes from this tree process. 

The second one is to dress the photon shower to the process Eq. flTBD in the similar 
way as the third case in the previous section, It means that the two hard photons 
are generated by QEDPS. 

The third one is to apply QEDPS to the process Eq. (|T9|) . In this case, two hard 
photons originate from the hard process Eq. and others from the simulation are 
not observed. 

The cuts for the hard photon is the same as before. Table 4 shows that even 
though the double-photon is required, QEDPS works in a consistent way indepen- 
dent of the allocation of the hard photon. 



V^(GeV) 






HZ-f-f/ps 


200 


(2.120±0.004)xl0~ 3 


(1.644±0.018)xl0~ 3 


(1.640±0.003)xl0~ 3 


300 


(2.161±0.004)xl0~ 3 


(2.040±0.017)xl0~ 3 


(1.922±0.004)xl0~ 3 


400 


(1.384±0.003)xl0- 3 


(1.370±0.011)xl0- 3 


(1.257±0.003)xl0- 3 


500 


(9.481±0.014)xl0~ 4 


(9.883±0.062)xl0- 4 


(8.707±0.022)xl0- 4 



Table 4 Total cross sections(pb) : double-hard photon required. 



Fig.4 shows the comparison in the distributions. In this figure, the histograms are 
for HZ/'yy (tree) and the cross and the circle correspond to HZ/'yyps and HZ'yy/ps, 
respectively. There appears a discrepancy in the distributions of Q and Q z between 
the Born and the other two with QEDPS because the higher order corrections are 
missing in the first case. On the other hand, in the region of high Q, HZ/'yyps 
and HZ'j'-f/ps are quite consistent. In the low region of Q, however, there seems a 
difference. Since there is no perturbative calculation with higher order we cannot say 
definitely which distribution is better but it is known in the study of QCD parton 
shower that the LL approximation always gives bigger estimation than the exact 
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one at larger momentum transfer region. E 1 and Ej distributions agree well among 
three cases. As a conclusion, these results show the self-consistency of the QEDPS 
model again even for the double hard photons. 



5 Summary and discussions 

In this work we have tested QEDPS, a generator for the radiative corrections in e + e~ 
annihilation. The total cross section corrected by QEDPS agrees with that calcu- 
lated by the analytic formula of the structure function with the accuracy less than 
around 0.3% in the energy range of LEP200 to linear colliders. We checked the hard 
photon distribution of QEDPS against the calculation of the matrix element with 
the 0(a) corrections and got an agreement. In order to make the self-consistency 
check, we applied QEDPS not only to e + e~ — > HZ but also to e + e" — > HZ ^ 
and ifZ°77. In these cases we have introduced the virtuality cut to avoid the 
double-counting of the hard photon. As the result it has been shown that QEDPS 
is self-consistent in various distributions. Through all of these studies we conclude 
that the model for the photon showers for the ISR is established fairly well even 
in the LL approximation. This is also demonstrated for more complicated process 



such as e + e~ — > yTv^ud as described in Ref. [13 . 

In the present work we have restricted the radiation only to ISR. The extension 
to include the final state radiation(FSR) is straightforward. In this case, however, 
a problem emerges how to incorporate with the interference between ISR and FSR, 
which is, in general, not so small to be ignored safely. Also the extension to the 
next-to-leading-logarithm(NLL) is urgent to make more precise predictions. This 
can be done in the same way as that in QCD parton shower These problems 
are now under investigation. 
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Figure Captions 

Fig.l Comparison of the distributions between H Z7 (tree) and HZ'-f with up to 
the 0(a) corrections at v / s=500GeV. The plots show the former and the histograms 
the latter. Q and Q z are the virtuality and the longitudinal momenta of the s-channel 
virtual boson, respectively. E 1 and E~{ are the energy and the transverse momentum 
of the required hard photon with respect to the beam axis, respectively. 

Fig. 2 Comparison of the distributions between HZ'-f with up to the 0(a) cor- 
rect ions (histograms) and iJZ/7ps(plots) at v /s=500GeV. 

Fig.3 Comparison of the distributions between HZ'-f with up to the 0(a) cor- 
rections(histograms) and HZ^ /ps(p\ots) at v / s=500GeV. 

Fig.4 Comparison of the distributions among HZ'yy (histogram), iJZ/77ps (star 
plots) and HZ'yy/ps (circle plots) at v /5=500GeV. The variables Q, Q z , E~' and E2 
are the same as in Fig. 1. 
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Figure 3: 
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